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SUPERNOVA FEEDBACK AND 
THE BEND OF THE TULLY-FISHER RELATION 
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RESUMEN 

Estudiamos el origen de la relation de Tully-Fisher analizando simulaciones hidrodinamicas en un universo 
ACDM. Encontramos que las galaxias mas pequenas exhiben masas estelares menores que aquellas predichas 
por el ajuste lineal de las galaxias mas masivas, en consistencia con las observaciones. En este modelo, estas 
tcndcncias son generadas por la action mas eficicnte de la retroalimcntacion al medio por Supernovas en la 
regulation de la formation estelar de las galaxias mas pequenas. Sin introducir parametros dependicntes de 
escala, el modelo prcdice un cambio de pendicnte en la relation de Tully-Fisher para una velocidad caracterfstica 
de ~ 100 km s , en acuerdo con resultados observacionales y tcoricos previos. 

ABSTRACT 

We have studied the origin of the Tully-Fisher relation by analysing hydrodynamical simulations in a ACDM 
universe. We found that smaller galaxies exhibit lower stellar masses than those predicted by the linear fit to 
high mass galaxies (fast rotators), consistently with observations. In this model, these trends are generated 
by the more efficient action of Supernova feedback in the regulation of the star formation in smaller galaxies. 
Without introducing scale-dependent parameters, the model predicts that the Tully-Fisher relation bends at a 
characteristic velocity of ~ 100 km s , in agreement with previous observational and theoretical findings. 
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1. INTRODUCTION 

The Tully-Fisher relation is very important for 
constraining galaxy formation models because it 
links two fundamental properties of galaxies such as 
the stellar (sTFR) or baryonic (bTFR) mass and the 
depth of the potential well. 

In particular, there are observational evidences 
that the sTFR does not follow a linear trend for 
the whole range of obs erved rotation velocities . Ac- 
cording to the results of lMcGaugh et al. 1 (|2000| ), the 
sTFR bends at ~ 90 km s _1 in such a way that 
smaller galaxies have lower stellar masses than those 
derived from the ex trapolation of the l inear fit to fast 
rotators. Moreover. McGaugh et al. ( 2010h have re- 
ported that there is also a bend in the bTFR but at 
a lower rotation velocity (~ 20 km s _1 ). 

2. SIMULATIONS 
In this work, we studied the role of Super- 
nova (SN) feedback on the shape of the sTFR and 



bTFR by analysing hydrodynamical simulations in 
a ACDM universe Q. A version of the chemical 
code GADGET-3 including treatments for metal- 
dependent radiative c ooling, stochastic star fq rma- 



tion and SN feedback (jScannapieco et al. 1120061 ) was 
employed. The simulated volume corresponds to a 
cubic box of a comoving 10 Mpc Yr 1 side length. The 
simulation has a mass resolution of 6 x lO 6 M /i _1 
and 9 x lOpM^hr 1 for the dark and gas phase, re- 
spectively. 

Simulated disc-like galaxies we re identified fol 
lowing the methods describe by iDe Rossi et al. 
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((20101). The mean properties of galactic systems 
were estimated at the baryonic radius i?b a r, defined 
as the one which encloses 83 per cent of the baryonic 
mass of the systems. We found that the tangential 
velocity of these systems constitutes a good repre- 
sentation of their potential well so that, for the sake 
of simplicity, we used the circular velocity estimated 
at i?bar as the kinematical indicator for our study. 

3. RESULTS AND DISCUSSION 

The massive-end of the simulated sTFR can be 
fitted by a power-law of the form log(M* / 'M©/i _1 ) = 
(3.68 ± 0.09) x log(V/ 100 km s" 1 ) +(9.42 ± 0.26), 
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in general good agreement with observations. How- 
ever, at rotation velocities below ~ 100 km s _1 , the 
sTFR becomes steep and the residuals of the linear 
fit depart syste matically from zero, con sistently with 
the findings of iMcGaugh et al~ ( 2000l ). To analyse 
the role of SN feedback on the origin of the bend 
of the sTFR, we run a feedback-free simulation. We 
found that, when SN feedback is suppressed from the 
model, the sTFR describes a linear behaviour indi- 
cating the crucial role of SNe in the modulation of 
this relation. 

With respect to the bTFR, both the SN feedback 
model and the feedback-free run show a single slope 
for the relation at least for the range of velocities 
resolved by these simulations (40 km s _1 < V < 
250 km s _1 ). Thes e results are not in disa greement 
with the findings of McGaugh et al. I ( 2010l ) because 
they reported a bend in the bTFR at ~ 20 km s _1 . 

By analysing the simulations at higher rcdshifts 
(0 < z < 3), we obtained similar trends. We have 
also checked that our results are robust against nu- 
merical res olution. For fur ther details, the reader is 
referred to iDe Rossi et al. 



(l2010h . 



Finally, we explore how the SN feedback model 
is capable of reproducing these behaviours without 
introducing scale-dependent parameters. To anal- 
yse the impact of galactic outflows in the simulated 
galaxies, we defined the fraction /™ as the ratio be- 
tween the baryonic mass within the virial radius to 
the one inferred from the universal baryonic fraction 
(fib/^m)- In these simulations, /™ is within the 
range 0.2 — 0.8 showing that for the whole sample 
a significant percentage of the gas in blown away as 
a consequence of efficient galactic winds. The more 
prominent losses are obtained for smaller galaxies. 

Moreover, by comparing the SN feedback model 
with the feedback-free run, we found that SN winds 
generate an important decrease in the star formation 
activity of galaxies with the larger effects in slow- 
rotating systems, consistently with the bend of the 
sTFR. Indeed, our results suggest that the role of 
SN feedback on the regulation of the star formation 
strongly depends on the mass of the galaxies. In 
particular, we distinguish two distinct regimes for 
the thermodynamical transitions of the gas phase. 
In the case of smaller galaxies, the virial temper- 
atures are lower and SN heating is more efficient 
at promoting gas from a cold to a hot phase (see 
IScannapieco et al~ 20061 for more details about the 
model). However, the cooling times of these systems 
are shorter than the dynamical times and the hot 
gas can return to the cold phase in short time-scales. 
Therefore, for slow-rotators, SN feedback leads to a 



self-regulated cycle of heating and cooling strongly 
influencing the star formation activity of these sys- 
tems. In the case of massive galaxies, the hot phase 
is established at a higher temperature and SN heat- 
ing cannot generate an efficient transition of the gas 
from the cold to the hot phase. Meanwhile, the 
cold gas remains available for star formation. On 
the other hand, the cooling times for these galaxies 
get longer compared to the dynamical times and the 
hot gas is able to remain in the hot phase during 
longer time-scales. Hence, SN feedback is not effi- 
cient at regulating the star formation in larger galax- 
ies. Interestingly, in this model, the transition from 
the efficient to the inefficient cooling regime for the 
hot-gas phase occurs at the same characteristic rota- 
tion velocity where the sTFR bends (~ 100 km s _1 ) 
and i s also in agreement with previous observational 



e.g. IMcGaugh et al. I|2000D and theoretical works 



(iLarson Ill974t iDekel fc Silk Ill986h . 



4. CONCLUSION 

We studied the Tully-Fisher relation by perform- 
ing numerical simulations in a cosmological frame- 
work. We obtained a sTFR and a bTFR in gen- 
eral agreement with observations. The SN feedback 
model is able to reproduce the observed bend of the 
sTFR at a characteristic velocity of ~ 100 km s™ 1 
without introducing scale dependent parameters. 
We found that this characteristic velocity separates 
two different regimes: slow-rotators develop a self- 
regulated cycle of efficient SN heating and radiative 
cooling while, in fast-rotators, the cold an hot gas- 
phases see m to be more disconnect ed. The reader is 
referred to De Rossi et al. ] (|2010h for more details 
about this work. 
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